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We consider the problem of diffusion-controlled evolution of the system A-particle island - B- 
particle island at propagation of the sharp annihilation front A + B — » 0. We show that this general 
problem, which includes as particular cases the sea-sea and the island-sea problems, demonstrates 
rich dynamical behavior from self-accelerating collapse of one of the islands to synchronous expo- 
nential relaxation of the both islands. We find a universal asymptotic regime of the sharp front 
propagation and reveal limits of its applicability for the cases of mean-field and fluctuation fronts. 
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For the last decades the reaction-diffusion system A + 
_B — > 0, where unlike species A and B diffuse and anni- 
hilate in a (i-dimensional medium, has acquired the sta- 
tus of one of the most popular objects of research. This 
attractively simple system, depending on the initial con- 
ditions and on the interpretation of A and B (chemical 
reagents, quasiparticles, topological defects, etc.), pro- 
vides a model for a broad spectrum of problems [l| , . 
A crucial feature of many such problems is the dynami- 
cal reaction front - a localized reaction zone which prop- 
agates between domains of unlike species. 

The simplest model of a reaction front, introduced al- 
most two decades ago by Gain and Racz (GR) 0], is a 
quasi-one-dimensional model for two initially separated 
reactants which are uniformly distributed on the left side 
(x < 0) and on the right side (x > 0) of the initial 
boundary. Taking the reaction rate in the mean-field 
form R(x,t) = ka(x,t)b(x,t), GR discovered that in the 
long-time limit kt — > oo the reaction profile R(x,t) ac- 
quires the universal scaling form 

where xj oc t 1 / 2 denotes the position of the reaction front 
center, Rf oc t~P is the height, and w oc t a is the width 
of the reaction zone. Subsequently, it has been shown 
[H]-[H that the mean field approximation can be adopted 
at d > d c = 2, whereas in ID systems fluctuations play 
the dominant role. Nevertheless, the scaling law (1) takes 
place at all dimensions with a = 1/6 at d > d c = 2 and 
a = l/4at<i=l, so that at any d the system demon- 
strates a remarkable property of the effective " dynami- 
cal repulsion" of A and B: on the diffusion length scale 
Li) oc t 1 / 2 the width of the reaction front asymptotically 
contracts unlimitedly: w/Lrj — > as t — > oo. Based 
on this property a general concept of the front dynamics, 
the quasistatic approxmation (QSA), has been developed 
0, Hi @]j 0| which consists in the assumption that for 
sufficiently long times the kinetics of the front is gov- 
erned by two characteristic time scales. One time scale 
tj = —(dlnJ/dt) -1 controls the rate of change in the 
diffusive current J = J a = \Jb\ of particles arriving at 
the reaction zone. The second time scale tf oc w 2 /D 



is the equilibration time of the reaction front. Assuming 
that tf/tj <C 1 from the QSA in the mean-field case with 
D A , B = D it follows that 0, @, 0] 

R f ~ J/w, w ~ (D 2 /Jk) 1/3 , (2) 

whereas in the ID case w acquires the fc-independcnt 
form w ~ {Djjfl 2 Q, [|. On the basis of the QSA a 
general description of spatiotemporal behavior of the sys- 
tem A + B —> has been obtained for arbitrary nonzero 
diffusiviti es I l0l| which was then generalized to anomalous 
diffusion pT| . diffusion in disordered systems [12|], dif- 
fusion in systems with inhomogeneous initial conditions 
(l3| . and to several more complex reactions. Following 
the simplest GR model [3( the main attention has been 
traditionally focused on the systems with A and B do- 
mains having an unlimited extension, i.e., with unlimited 
number of A's and B's particles, where asymptotically 
the stage of monotonous quasistatic front propagation is 
always reached: tf/tj — > as t — > oo. 

Recently, in the work [lij a new line in the study of 
the A + B — * dynamics has been developed under the 
assumption that the particle number of one of the species 
is finite, i.e., an A particle island is surrounded by the 
uniform sea of particles B. It has been established that 
at sufficiently large initial number of A particles A^o and a 
sufficiently high reaction rate constant k the death of the 
majority of island particles N(t) proceeds in the universal 
scaling regime N = N g{t/t c ), where t c oc N§ is the 
lifetime of the island in the limit k, Nq — ► oo. It has been 
shown that while dying, the island first expands to a 
certain maximal amplitude x^ oc No and then begins to 
contract by the law Xf = x¥ £f(t/t c ) so that on reaching 
xy (the turning point of the front) 

t M /t c = l/e, N M /No = 0.19886... (3) 

and, therefore, irrespective of the initial particle number 
and dimensionality of the system w 4/5 of the particles 
die at the stage of the island expansion and the remaining 
~ 1/5 at the stage of its subsequent contraction. 

In this Rapid Communication we consider a much more 
general problem of the A + B — > annihilation dynamics 



2 



with the initially separated reactants under the assump- 
tion that the particle number of the both species is finite. 
More precisely, we consider the problem on diffusion- 
controlled death of ^4-particle and B-particle islands at 
propagation of the sharp annihilation front A + B — > 0. 
We show that this island-island (II) problem, of which 
particular cases are the GR sea-sea (SS) problem and 
the island-sea (IS) problem [l4|], exhibits rich dynamical 
behavior and we reveal its most essential features. 

Let in the interval x e [0, L] particles A with concen- 
tration a and particles B with concentration b be ini- 
tially uniformly distributed in the islands x € [0, £) and 
x G (i,L], respectively. Particles A and B diffuse with 
diffusion constants Da and Db and when meeting they 
annihilate A + B — ► with a reaction constant k. We 
will assume, as usually, that concentrations a(x, t), b(x, t) 
change only in one direction (flat front) and we will con- 
sider that the boundaries x = 0, L are impenetrable. 
Thus, our effectively one dimensional problem is reduced 
to the solution of the problem 

da/dt = D A V 2 a~ R, db/dt = D B V 2 b - R (4) 

in the interval x £ [0, L] at the initial conditions a(x, 0) = 
ao9(£ — x), b(x, 0) = bo&(x — I) and the boundary con- 
ditions V(a, b) \ x =o,l= where 9(x) is the Heaviside step 
function. To simplify the problem essentially we will as- 
sume Da = Db = D. Then, by measuring the length, 
time, and concentration in units of L, L 2 /D, and &0j re- 
spectively, i.e. assuming L = D = bo = 1, and defining 
the ratio of initial concentrations ao/&o = T and the ra- 
tio £/L — q, we come from (4) to the simple diffusion 
equation for the difference concentration s = a — b 



ds/dt = V 2 s, 



(5) 



in the interval x € [0, 1] at the initial conditions 

s (xe[0,q))=r, s (xe(q, !]) = -!, (6) 



with the boundary conditions 

Vs L=o,i= 0. 



(7) 



According to the QSA for large k — > oo at times t oc 
fc _1 — > there forms a sharp reaction front w/xf — > 
so that the solution s(x,t) defines the law of its propa- 
gation s(xf,t) — and the evolution of particle distri- 
butions a = s(x < Xf) and b = \s\(x > Xf). In the 
limits sea-sea [3|(£ — > oo,L — > oo) or island-sea [14[ [t 
finite, L — > oo) the corresponding solutions sss(x,t) and 
sis{x,t) describe the initial stages of the system's evolu- 
tion at times \fi <C g, 1 — q and gCv'fCl, respectively. 
The general solution to Eqs.(5)-(7) for arbitrary r, q and 
t has the form 



s(x, t) = A + ^2 A n( r , q) cos ( 



mix e 



(8) 



n=l 



where coefficients A n (r,q) = 2(r + 1) sin(n7rg)/n7r and 
A(r, q) — Na — Nb = rq — (1 — q) is the difference of 



the reduced number of A and B particles which remains 
constant. At t > l/ir 2 we find 



A + Ai(r, q) cos(wx)e 



(9) 



Taking s(xf,t) = we obtain from (9) the law of the 
front motion 

cos{TTx f ) = Ce 7r2 * + -- - , (10) 

where coefficient C can be represented in the form 

C = -A/A! = q{r, - r)/A 1 = (g, - q)/q*A 1 , (11) 

where q+ = \/(r + 1) and r* = (1 — q)/q are the critical 
values of q,r at which C reverses its sign. From Eq.(10) 
it follows that at \C\ < 1/e and r ^ r*, q ^ q+, when the 
ratio of the initial particle numbers 



= N A g = r_ = (1 - q*)q ^ ^ 
Nbo (1 - q)q* 



(12) 



the front x / (t) moves either towards the boundary x = 
(p < 1) or towards the boundary x = 1 (p > 1) so that in 
the limit k — > oo the island of a smaller particle number 
(A or B, respectively) dies within a finite time 



< c = (lA 2 )|ln|C||. 
From Eqs. (10) and (13) we obtain 

x f = (l/7r)arccos(±e 7r2( *- t '= ) ), 



(13) 



(14) 



(here and in what follows the upper sign corresponds to 
p < 1 and the lower sign corresponds to p > 1) whence 
for the front velocity Vf — if we find 



v f = -ircot(irxf) = iF7r/(Ve 27r2 (*<="*) - 1). (15) 

Making use then (13), for the distribution of particles 
(a = s(x < Xf),b = \s\(x > Xf) [TH) at p ^ 1 we obtain 



s = A(l q= cos(7ra)e 



7r 2 (tc-t)\ 



(16) 



Thus from the condition Na = J Q Xf sdx = Nb + A we 
find the laws of decay of the A and B particle number 



N A = (|A|/7r)(Ve 27r2 (*<=-*) -l^irx f ) (17) 

and then we derive finally the diffusive boundary current 
in the vicinity of the front 



J 



-ds/dx | B=a/ = 7r|A|Ve 2 - 2 ^-*) -1, 



(18) 



which according to (2) defines evolution of the amplitude 
Rf(t) and of the width of the front w(t). 

From Eqs. (13)-(18) we immediately come to the fol- 
lowing important conclusions: for arbitrary r and q at 
p < 1 or p > 1 (i) the motion of the front is the universal 
function of the " distance" to the collapse time t c — t with 
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the remarkable property Xf(t c — t) = 1— Xf(t c — t). More- 
over, the front velocity Vf is the unique function of Xf 
with the remarkable symmetry Xf <-> 1 — Xf,Vf <-> —Vf, 
(ii) the reduced particle number Na/\A\ and the re- 
duced boundary current J/|A| are universal functions 
of t c — t with the remarkable properties N^(t c — t) = 
N>(t c - 1) - |A| and J<(t c - t) = J>(t c - t). 

Introducing the relative time T = t c —t, from Eqs. (13)- 
(18) in the vicinity T <C 1/ir 2 of the critical point t c 
we come to the universal power laws of self-accelerating 
collapse (\v f \ oc T" 1 / 2 ) 

xf,l-x> = V2T+--- , (19) 
7V<,7V> = (\/8/3)7r 2 |A|T 3 / 2 + --- , (20) 

J = V2vr 2 |A|VT+--- . (21) 

At large i c ^> 1 /ir 2 far from the critical point T > I/71- 2 
according to Eqs. (13)-(18) there is realized the interme- 
diate exponential relaxation regime (\vj\ oc e _Tr r ) 

x f> = 1/2 t e - *" 7 /* + ■ • • , (22) 
7V<<> = (|A|/7rK 2r (l T ^ T /2 + •••), (23) 

J = ^|A|e" 2r (l-e- 2 " 2T /2 + ---), (24) 

which in the limit t c — > 00 (|C|, \ g — 1| — ► 0) becomes 
dominant. Thus, at large £ c 3> l/7r 2 the point a;/ f=a 
1/2 (stationary front) is an "attractor" of trajectories. 
Exactly at the critical point = 1 from Eqs. (9) and 
(10) we find x*j = 1/2 and obtain 

"■<»°-(i)w^r" + -< <25 » 

J+ = 2(sin(vrg)/g)e" 7r2t + • • • . (26) 

In order to answer the question of when and how the 
"attractor" = 1/2 is reached it is necessary to retain 
the next term (n = 2) in the sum (8). With allowance 
for the first two terms one can easily obtain 

x} = 1/2- D(q)e^ 2t + ■■■ (27) 

where D(q) = (A^/irAi) = wa.(2%q)/2Trsia.(irq). Ac- 
cording to (27) at q = 1/2 the coefficient D reverses 
its sign, therefore, as it is to be expected, at q < 1/2 
and q > 1/2 the front reaches the attractor x*f = 1/2 
from the left and the right, respectively. By combin- 
ing (22) and (27), at small but finite C we have xf = 



l/2-Ce" t/ir-Der^ *H . We thus conclude that un- 
der the condition DC > there arises the turning point 
of the front (vj 1 = 0) with the coordinates 

t M = (1/4tt 2 ) \n(X M \D/C\) + ■■■ , (28) 

xf = 1/2 - m M D\C/D\ 3 / 4 + ■ • • , (29) 

where Am = 3n,mM = 4/(37r) 3 / 4 , whereas at DC < 
there arises the inflection point of the front trajectory 
(|u?| = min|uf|) with the coordinates t s ,x S f which are 
determined by Eqs. (28), (29) with the coefficients 
A s = 3\M,m s = 2mM/ / (3) 3/ ' 4 . The analysis presented 
demonstrates the key points of evolution of the island- 
island system at arbitrary q and r with C(q,r) < 1/e. 
Below we will focus on a detailed illustration of this evo- 
lution from the initial island-sea configuration (q <C 1). 

A remarkable property of the island-sea configuration 
q <C 1 is that at r 3> 1 the A(p) = p — 1 value and all 
the coefficients A n (p) = 2p up to n oc 1/q ^> 1 become 
unique functions of p. Therefore, the system's evolution 
at t 3> q 2 is determined by the sole parameter p. At 
q 2 <C t <C 1 we have the scaling IS regime [l4| 

x f = %/27ln 1/2 (p 2 /7rf), t c (p) = p 2 /tt (30) 

with xf = p^2/ne, t M = p 2 jite. For t > l/4ir 2 with 
allowance for two principal modes (n = 1, 2) we obtain 
from (8) 

x f = (1/tt) arccos[G(p, t)e^ H j\\ (31) 

where G(p, t) = yjl + 8C{p)e~ 2 * 2t + 8 e - 67r2 * - 1 and 
C(p) = (1 — p)/2p. For the time of collapse t c (p) we 
derive from (8) the general equation for arbitrary p 

00 

^(±l)" e -" 2 - 2 Mp) =±\C(p)\ (32) 

n=l 

whence for the leading terms in accord with (31) we find 

t c (p) = (\\n\C\\±\C\ 3 + ---)/7r 2 . (33) 

Using small t representations of the series (32), one can 
easily show that, with the growing p, t c initially grows 
by the law t c (p) = p 2 (l + 4e~ 7r / p2 + •■■)/%, then it 
passes through the critical point t c (p*) — ► 00 accord- 
ing to Eq.(33), and finally at large p decays by the law 
t c (p) oc 1/ln/j. From Eqs. (31) and (17) for the starting 
points tM, s of front self-acceleration at small \C\ we find 

t M , s /t c = l/A + p MiS /\ln\C\\+--- (34) 

with the number of A particles N^' s /Nao °c | C | 1 / 4 where 
Pm = /3 s /2 = In 3/4. Remarkably, the same as for the 
scaling IS regime (3), (30) in the vicinity \p — pJ <C 
1 the ratio t^/tc reaches the universal limit tuj^c — 
1/4. In Fig. 1 are shown the calculated from (30) and 
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(31) trajectories of the front Xf(t), which illustrate the 
evolution of the front motion with the growing p. It is 
seen that to p ss 0.7 the death of the island A proceeds in 
the scaling IS regime (30) (t m jt c = 1/e) , then the x / (i) 
trajectory begins to deform, and at small \p — p±\ -C 1 
the regime of the dominant exponential relaxation (22)- 
(24), (34) (t M /t c ~ 1/4) is reached. After the critical 
point p^ = 1 has been crossed, the death of the island A 
is superseded by the death of the island B, so the front 
trajectory becomes monotonous, and the stopping point 



of the front 



0) "transforms" to the point 
(Vf = minuy cx 



of maximal deceleration of the front x s j , t 
|C| 3 / 4 ) which at large p shifts by the law 1 — x S f cx 1/ 
with t s cx 1/lnp. 

One of the key features of the island-island problem is 
a rapid growth of the front width w while the islands are 
dying. Therefore, to complete the analysis we have to 
reveal the applicability limits for the sharp front approx- 
imation i] = w/min(xf, 1 — Xt) <C 1. By substituting 
(21) into (2) we obtain for the self-accelerating collapse 
i] ~ (Tq/T)^ where for the mean- field front /kmf = 2/3 
and Tq 1F = 1/ ^/|A|fc. For perfect diffusion-controlled 
3D reaction k ~ Dr a where r a is the annihilation radius. 
Thus, as our k is measured in the units of D/L 2 b [3| 
for the dimensionless k we have k = r a L 2 bo. Substituting 
here r a ~ 10 -8 cm, L ~ 10 cm and bo ~ 10 22 cm~ 3 we 
find k ~ 10 16 and derive T^ F ~ lO^/VlA] so that for 
not too small |A| (|p — p*| ^s> 10 -8 ) the sharp front is not 
destroyed almost down to the point of collapse. Clearly 
that at small |A| — ► the "destruction" of the front has 
to occur already at the stage of exponential relaxation 



(22)-(26). Substituting (26) into (2) for the exponen- 
tial relaxation we find r\ ~ e 1 " 1 (* _t «) where ^mf = 1/3 
and tW = (lnfc)/7r 2 . Substituting here k ~ 10 16 we ob- 



MF 



3.7 and then from (25) we find Nf* (77 



tain fc 

0.1) /Nq ~ 10^ 13 . The analogous calculation for the fluc- 
tuation ID front gives p F = 3/4, - l/(|A|n ) 2/3 and 
v F = 1/2, tq — (lnno)/7r 2 where no = Lb$. Substitut- 
ing here n - 10 6 we find Tj - 10~ 4 /| A| 2 / 3 , t% - 1.4 
and n F (r] = 0.1)/no ~ 10~ 4 . We conclude that both for 
the MF and the fluctuation fronts the vast majority of 
the particles die in the sharp front regime, therefore, the 
presented theory has a wide applicability scope. 

In summary, the evolution of the system island A- 
island B at the sharp annihilation front A + B — ► prop- 
agation has been first considered and a rich dynamical 
picture of its behavior has been revealed. The presented 
theory may have a broad spectrum of applications, e.g. 
in description of electron-hole luminescence in quantum 
wells 15 1, formation of nontrivial Liesegang patterns [lj|, 
and so on. Of special interest is the analogy of the island- 
island problem with the problem of annihilation on the 
catalytic surface of a restricted medium where for un- 
equal species diffusivities in a recent series of papers (l7| 
the phenomenon of annihilation catastrophe has been dis- 
covered. Study of the much more complicated case of 
unequal diffusivities and comparison with the annihila- 
tion dynamics on the catalytic surface is a generic and 
challenge problem for future. 
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FIG. 1: (Color online) Evolution of the front trajectories 
Xf(t) with growing p, calculated from Eqs. (30) (blue lines) 
and (31)(red circles) at p = 0.5,0.7,0.9,0.98, 1, 1.02, 1.1 and 
2 (from left to right). The region of the scaling IS regime is 
shaded. 
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